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Abstract 

The focusing of the radiation generated by a polarization current with a superluminally rotating 



O ' distribution pattern is of a higher order in the plane of rotation than in other directions. Conse- 

c/3 , quently, our previously published asymptotic approximation to the value of this field outside the 

equatorial plane breaks down as the line of sight approaches a direction normal to the rotation 
^ ■ axis, i.e., is nonuniform with respect to the polar angle. Here we employ an alternative asymptotic 

^ : 

. expansion to show that, though having a rate of decay with frequency (fi) that is by a factor of 

m 

I order fi'^^'^ slower, the equatorial radiation field has the same dependence on distance as the non- 
spherically decaying component of the generated field in other directions: it, too, diminishes as the 

. inverse square root of the distance from its source. We also briefly discuss the relevance of these 

- ^ , results to the giant pulses received from pulsars: the focused, nonspherically decaying pulses that 

X 

' arise from a superluminal polarization current in a highly magnetized plasma have a power-law 

a ! 

spectrum {i.e., a flux density S oc whose index (a) is given by one of the values —2/3, —2, 
-8/3, or -4. 
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I. INTRODUCTION 



Radiation by polarization currents whose distribution patterns move faster than hght 
in vacuo has been the subject of several theoretical and experimental studies in recent 
years l|-l8|. When the motion of its source is accelerated, this radiation exhibits features 
that are not shared by any other known emission. In particular, the radiation from a 
rotating superluminal source consists, in certain directions, of a collection of subbeams 
whose azimuthal and polar widths narrow (as Rp^ and -Rp^, respectively) with distance Rp 
from the source [sl. Being composed of tightly focused wave packets that are constantly 
dispersed and reconstructed out of other waves, these subbeams neither diffract nor decay 
in the same way as conventional radiation beams. The field strength within each subbeam 
diminishes as Rp^^"^, instead of -Rp\ with increasing Rp jo-^. 

In earlier treatments 0, Q] , we evaluated the field of a superluminally rotating extended 
source by superposing the fields of its constituent volume elements, i.e., by convolving its 
density with the familiar Lienard-Wiechert field of a rotating point source. This Lienard- 
Wiechert field is described by an expression essentially identical to that which is encountered 
in the analysis of synchrotron radiation, except that its value at any given observation time 
receives contributions from more than one retarded time. The multivalued nature of the 
retarded time gives rise to the formation of caustics. The wave fronts emitted by each 
constituent volume element of a superluminally moving accelerated source possess a cusped 
envelope on which the field is infinitely strong (see Figs. 1 and 4 of Ref. Q]). Correspondingly, 
the Green's function for the problem is nonintegrably singular for those source elements 
that approach the observer along the radiation direction with the speed of light and zero 
acceleration at the retarded time (see Fig. 3 of Ref. j^): the cusp of the envelope of wave 
fronts emanating from each such element is a spiraling curve extending into the far zone that 
passes through the position of the observer. When the source oscillates at the same time 
as rotating, the Hadamard finite part of the divergent integral that results from convolving 
the Green's function with the source density has a rapidly oscillating kernel for a far-field 
observation point. The stationary points of the phase of this kernel turn out to have different 
orders depending on whether the observer is located in or out of the equatorial plane. 

To reduce the complications posed by the higher-order stationary points of this phase, 
we restricted the asymptotic evaluation of the radiation integral thus obtained in Refs. 3, Q] 
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to observation points outside the plane of rotation, i.e., to spherical polar angles 6p that do 
not equal 7r/2. The purpose of this paper is to evaluate the field of a superluminally rotating 
extended source also for the smaller class of observers at polar coordinate dp = 7i/2 and 
to obtain, thereby, a more global description of the nonspherically decaying radiation that 
is generated by such a source. The asymptotic expansion presented in Refs. 0, S] breaks 
down in the case of a subbeam that is perpendicular to the rotation axis because there is 
a higher-order focusing associated with the waves emitted by those source elements whose 
actual speeds (rather than the line-of-sight components of their speeds) equal the speed of 
light as they approach the observer with zero acceleration. 

Here, we present a brief account of the background material on the radiation field of a 
rotating superluminal source in Section 2, and the asymptotic evaluation of this field for an 
equatorial observer in Section 3. In Section 4, we give a description of the spectral properties 
of the nonspherically decayingcomponent of this radiation in the light of the present results 
and those obtained in Refs. 7, ^, and discuss the relevance of these properties to pulsar 
observations. 



II. BACKGROUND: RADIATION FIELD OF A ROTATING SUPERLUMINAL 
SOURCE 

We base our analysis on a generic superluminal source distribution 7|, ISj, which has been 
created in the laboratory [2;]. This source comprises a polarization current density j = dP/dt 
for which 

Pr,<p,z{r, (f, Z, t) = Sr^Lp^zij^ z) COs{nup) COs(fit), — TT < < TT, (l) 

with 

(f = If — ut, (2) 

where Pr,ip,z are the components of the polarization P in a cylindrical coordinate system 
based on the axis of rotation, s(r, z) is an arbitrary vector that vanishes outside a finite 
region of the (r, z) space, and m is a positive integer. For fixed t, the azimuthal dependence 
of the density ([T]) along each circle of radius r within the source is the same as that of a 
sinusoidal wave train, of wavelength 27rr/m, whose m cycles fit around the circumference 
of the circle smoothly. As time elapses, this wave train both propagates around each circle 
with the velocity ru and oscillates in its amplitude with the frequency Q. This is a generic 



source: one can construct any distribution with a uniformly rotating pattern Pr^^p^zif, 0, z) 
by the superposition over m of terms of the form Sr,^,z{f^, z, m) cos{m(p). 
The electromagnetic fields 

E = -VpA°-9A/(9(ctp), B = VpxA (3) 

arising from such a source are given, in the absence of boundaries, by the following classical 
expression for the retarded four-potential: 

^'^(xp, tp) = c-^ j d?x&tf{^, t)S{tp - t - R/c)/R, /i = 0, ■ ■ ■ , 3. (4) 

Here, (xp,tp) = {rp,{pp, zp,tp) and (x, t) = {r,{p,z,t) are the space-time coordinates of 
the observation point and the source points, respectively, R stands for the magnitude of 
R = xp — X, and fi = 1,2,3 designate the spatial components, A and j, of A'^ and in a 
Cartesian coordinate system. 

To find the retarded field that follows from Eq. (j4]) for the source described in Eq. ([1]), 
we first calculated in Ref. 7|] the Lienard-Wiechert field arising from a circularly moving 
point source with a speed ru > c, i.e., a generalization of the synchrotron radiation to 
the superluminal regime. We then evaluated the integral representing the retarded field 
of the extended source ([T]) by superposing the fields generated by the constituent volume 
elements of this source, i.e., by using the generalization of the synchrotron field as the 
Green's function for the problem. In the superluminal regime, this Green's function has 
extended singularities arising from the constructive intereference of the emitted waves on 
the envelope of wave fronts and its cusp. 

Labeling each element of the extended source ([T]) by its Lagrangian coordinate if and 
performing the integration with respect to t and (or equivalently ip and 0) in the multiple 
integral implied by Eqs. we showed in Ref. 7| that the resulting expression for 

the radiation field B (or E) consists of two parts: one that decays spherically (as Rp^, 
as in a conventional radiation field) and another, B'^^ (with E"*^ = fixB"*^), that decays 
nonspherically (as -Rp ) within the conical shell arcsin(l/f>) < ^p < arcsin(l/f<;) in the 
far zone. Here, {Rp,6p,ipp) are the spherical polar coordinates of the observation point P, 
f stands for ruj/c, n = R/i? is a unit vector in the radiation direction, and f< > 1 and 
f> > f< denote the radial boundaries of the support of the source density s. 

The expression for the nonspherically decaying component of the field within this conical 



shell, in the far zone, is 



(5) 



B- ^ - |iexp[i(r]/^)(^p + 37r/2)] Z,=,^ /iexp(-i/x^p) 
x\~^gj / fdf A~"'^/^Uj exp(— i/i0_), 

where /i± = (fi/cj) ± m, i^p = (/9p — utp, 

Qj = {1 —ifl/uj iQ/u), (6) 

ui = cos6'pe|| + s<^ej_, U2 = — cos6'pe|| + Sre_L, U3 = — sin 6'pe||, (7) 
A = (r^-l)(f2-l)-(i-ip)2, (8) 

(l)± = R± + (p±- (pp, (9) 
= (^p + 27r - arccos[(l T A^/^)/{ffp)], (10) 



and 



= p - ip)2 + f2 + f2 - 2(1 T AV2)]1/2; (11) 



see Eq. (47) of Ref. p]. In this expression, (f, z; fp, ;zp) stand for {ru/c, zu/c; rpu/c, zpuj/c), 
and e|| = e^xn/lexXfil (which is parallel to the plane of rotation) and e_L = nxey comprise a 
pair of unit vectors normal to the radiation direction fi (e^ is the base vector associated with 
the coordinate z) . The domain of integration consists of the part of the source distribution 
s(r, z) that falls within A > (see Fig. 4 of Ref. fl]). 

Both derivatives, d(f)^/df and d(f)^/dz, of the function 0_(f, z) that appears in the phase 
of the integrand in Eq. ([5]) vanish at the point f = 1, z = Zp, where the cusp curve of the 
bifurcation surface is tangent to the light cylinder (see Figs. 3 and 4 of Ref. |8j). However, 
d'^cj)-/ dr'^ diverges at this point, so that neither the phase nor the amplitude of the kernel 
of the integral in Eq. ([5]) are analytic at f = 1, z = 5p. Only for an observer who is located 
outside the plane of rotation, i.e., whose coordinate zp does not match the coordinate z of 
any source element, is the function 0_(r, z) analytic throughout the domain of integration. 
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To take advantage of the simplifications offered by the analyticity of 0_ as a function of r, 
we restricted the analyses in Refs. 7|, |8| to observation points for which 6p ^ 7r/2. 
In the calculation that follows, we find an asymptotic approximation to the integral 

T= / fdfdi A"^/^Ujexp(-i/i0_) (12) 
Ja>o 

in Eq. ([5]) that is valid in the plane of rotation, i.e., for 6p = 7r/2. We shall treat only the 
case of positive fi; T{ii) for negative can then be obtained via X(— /i) = X(/i)*. 

III. ASYMPTOTIC VALUE OF THE FIELD FOR AN EQUATORIAL OBSERVER 
IN THE FAR ZONE 

Since the main contribution toward the value of the field at = 7r/2 is made by the 
source elements that lie in the vicinity of the critical point f = 1, i = zp, the first step in the 
asymptotic evaluation of B"^^ is to replace (f, z) by a new pair of variables (p, a) for which 
the phase function (j)_{p,a) is rendered analytic at this point: 

f = (l + p'coshV)l/^ (13) 

2 = zp + {rl- 1)^/Vsinhcr (14) 
This transformation replaces fA'^^^drd^ by pcoshcrdpdo" and yields 

0_(p, a) =[fl - 1 - 2{fl - Xfl'^p + {jI sinh^ a + l)pT^^ + Svr 
- arccos{fp^(l + p^ cosh^ (yY^'\^ + (^p - 



(15) 



which is analytic at p = a = 0. In the plane of rotation, i.e., for cr = 0, the two critical 
points designated as C and S in Refs. 7, ^ coalesce, and both derivatives, dcjy^/dp and 
(9^0_/9p^, of the resulting function 0_(p, 0) vanish at p = [see Eqs. (fTSl) and (!24l) below], 
so that the function (f)-{p, cr) is stationary at p = a = 0. 

To see that applying the method of stationary phase to the integral in Eq. ([5]) results in 
a valid asymptotic approximation for large Rp, let us begin by casting the a dependence 
of the phase (/)_ into a canonical form 3]. Since a = is an isolated stationary point of 
0_ (when regarded as a function of the single variable a), we may employ the following 
transformation: 

0_=0_|^=o + |&C', (16) 



in which 

p^[rl - 1 - (f|. - l)V2p + f^p^^ 



b 



da' 



a=0 (l + p2)[(f2 -1)1/2 _p] • (17) 

Equation f|T6l) expresses a as a function of C imphcitly. Repeated differentiations of this 
equation with respect to ( result in 

K. (18) 



da dC 



and so on, which when evaluated at C = supply the coefficients da/d(\u=o, d'^a/d(^\u=o, 
etc., in the Taylor expansion of a in powers of (. 

The integral X in Eq. (fT2|) can therefore be written as 

1 = j dpdCg(p,C)exp(-i/3C'), (20) 

where 

Q{pX) = pcoshcrUjexp(-i/x0_|^=o)c^cr/<9C, (21) 

da 6C^-(1 + p^coshV) ^ ..2 1 u/2 ,-22 u2 i-i /oo^ 
P^sinhacosha [^p " ^ " (^p " 1) ^ P + W cosh , (22) 

and /5 = The limits of integration are determined by the image of A > under 

transformation (fT6!) . 

The parameter h that multiplies the phase of the integrand in Eq. ( l20l) has a large value 
in the far zone: 

h - p'^fp, Rp > 1, (23) 

[see Eq. ( ITTl) ]. The asymptotic value of the integral X for large -Rp therefore receives its 
leading contribution from the immediate vicinity of C = 0, where the phase of its integrand 
is stationary [ot]. Replacing Q{p,C) Eq. ( l20l) by Q{p,0) and extending the range of 
integration with respect to ( to (—00, 00), we obtain 

X ~ {271 / p)'/^ j dppu,U=ori/2exp[-i(;U0_|,=o + n/A)]{da/dC)a=o 

~ (27r/p)i/2fpi/2 exp{-i[p(rp + 37r/2) + 7r/4]} (24) 
X y dpuj|^^(i_^p2)i/2^£=£p exp[— ip(arctanp — p)], .Rp 3> 1, 



where the integation extends over all values of p for which the source density s|j:^(]^^p2)i/2^^=£p 
is nonzero [see Eq. ([7j)]. In deriving this expression, we have inferred the value da /dC,\a=f) = 1 
of the indeterminate Jacobian that appears in (^(p, 0) from Eq. fll9p [or, equivalently, from 
Eq. (1221) and THopital's rule], and expressed (j)_\^^Q and h in terms of their far-field values 
by means of Eqs. (IT^ and 

The contribution B"^*^ toward the magnetic field B of the radiation is made by those volume 
elements of the source that approach the observation point P, along the radiation direction, 
with the speed of light and zero acceleration at the retarded time, i. e., by the source elements 
for which A = 0. Hence, the amplitude of the integrand in Eq. ([5]) has already been 
approximated by its leading term in powers of A^/^ = (fp — l)^^^p (see Refs. j?, Q]). To be 
consistent, we must also approximate the amplitude of the integrand in Eq. fl2^ by its value 
for p <^ 1: 

I ~(27r/p)i/2rp'/'u,-|,=i,,=,p exp{-i[p(fp + 3n/2) + 7r/4]} 

Hf|-i)i/2 (25) 
X / dp exp[— ip(arctanp — p)]. 



Jo 

where f> denotes the radial extent of the support of the source density s. This reduces to 

J = 3-'/'T{l){27r'f'fi-'/%\r=i,,=,, exp{-i[p(fp + 37r/2) + 7r/12]}f-^/' (26) 

in the regime p ^ 1, where we can approximate arctanp — p in the argument of the 
exponential by — and replace the upper limit of integration by oo [9|. 

Thus, Eqs. ([5]), ( |T2l) and ( l25l) jointly yield the following expression for the leading term 
in the asymptotic expansion, for Rp ^ 1, of the magnetic field of the radiation close to the 
plane 9p = n/2: 

B ~ - -i(27r)i/2^p^/^ csc^/^ exp[i(fi/cj)((/?p + 37r/2)] 
3 



(27) 



where 



X l/^r^^sgn(p)exp{-i[p((^p + fp + 37r/2) + ^sgn(p)]} 

3 

Hf 2 -1)1/2 

J= I dp exp[— ip(arctanp — p)]; (28) 



for, the contribution B"^*^ toward the magnetic field B of the radiation is larger by a factor 

^1/2 

of the order of R-p than the spherically decaying contribution. This is the counterpart of 
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Eq. (55) of Ref. [3] and Eq. (61) of Ref.Q] (the electric-field vector of this radiation is given 
by fixB as in any other radiation). 

Note that the remaining integral in the above expression reduces to 

J ~ 3-2/3^(1) exp(i7r/6)/i-^/=' (29) 
in the limit l/x] > 1 [see Eq. fl26D ]. 



IV. SPECTRUM OF THE NONSPHERICALLY DECAYING RADIATION: REL- 
EVANCE TO PULSAR OBSERVATIONS 

Eq. fl7r|) shows that the radiation field of a rotating superluminal source diminishes as 
Rp with the distance Rp also in the equatorial plane 9p = 7r/2. This differs from the 
corresponding result for 6'p 7^ 7r/2 [Eq. (55) of Ref. 7] mainly in its dependence on frequency. 
The Fourier transform s in Eq. (57) of Ref. Tj has the asymptotic dependence fi~^ on fi for 
a source density s(r, z) that is a smooth function of z. Therefore, when s{z) is smooth and 
the radiation frequency \fiuj\ appreciably exceeds the rotation frequency u, the field in the 
plane of rotation decays more slowly with frequency, by a factor of order /x^'^^, than does the 
field outside this plane. 

Since the azimuthal width of the generated subbeams (and hence the duration of the 
narrow signals that constitute the overall pulse) is independent of frequency, the flux density 
S of such signals {i.e., the power propagating across a unit area per unit frequency) is 
proportional to |Bp/AyU, where A/i ~ is the bandwidth of the radiation. The flux density 
of the emission described by Eq. f l271) thus depends on frequency as S* oc /i~^/^g||s(/i)p for 

^ 1. Here, |s(yu)| designates the frequency dependence of the factor s that enters the 
expression for the polarization P and the definitions of [see Eqs. ([1]) and ([7])]. The flux 
density of the corresponding emission outside the equatorial plane depends on frequency 
as S* oc /i~^g||s(yu)p for ^ 1 since, apart from the dependence s(/x) of a mutiplicative 
factor (such as electric susceptibility) in s, the Fourier transform s would decay as /x^^ for 
a smooth s{z). 

To compare the predictions of Eq. (p71) [and Eq. (55) of Ref. 7] with the observed spectra 
of the giant pulses from pulsars, we therefore need to estimate the frequency dependence 
of the electric susceptibility (contained in the factor s) for the magnetospheric plasma of 
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a pulsar. The simple classical model of propagation of electromagnetic disturbances in a 
cold magnetized plasma yields a dielectric tensor, and hence an electric susceptibility, whose 
components decay with frequency as (/icu)""^ when the frequency /icu of the disturbance that 
polarizes the medium is much lower than the gyration frequency of the electrons in the 
magnetized plasma; see, e.g., Eq. (7.67) of [l^. For a magnetic field as strong as that of a 
pulsar (~ 10^^ G), the Larmor frequency of an electron exceeds the highest radio frequencies 
at which the pulses are observed by a factor of order 10^, so that s(/i) oc /i^^ for pulsars. 

Using this result, we obtain S oc (fjii~^/^ for O-p = n/2, and S oc qj/J^''^ for 9p ^ 7r/2. 
Depending on whether the modulation frequency Q in the expression for qj [Eq. (Q] is com- 
parable to or much smaller than the frequency mu of the sinusoidal wave train characterizing 
the spatial distribution of the source, therefore, the spectral density of the nonspherically- 
decaying radiation is given by 

S (X fi 

-2/3^ ep = n/2, n/uj-\fi\, (30) 



Socfx-^, 9p^n/2, n/uj-\fx\, (31) 



S oc Op = 7i/2; Q/u < or j = 1, (32) 



or 



S oc ^ 7r/2; n/u < or j = 1. (33) 

In other words, the spectral index of the pulses portraying the subbeams can have any of 
the values -2/3, -2, -8/3, or -4. 

The range of spectral indices (—4 < a < —2/3) implied by Eq. fl27|) and its counterpart, 
Eq. (57) of Ref. j?], is consistent with that which characterizes the observed power-law 



spectra of the giant pulses from pulsars |lll4l3l| . For radio pulsars, the rotation frequency uj 



of the distribution pattern of the radiating polarization current is of the order of 1 rad/sec, 
and the oscillation frequency fiuj/2TT of the source density, of the order of 100 MHz, so that 
/i has a large value of the order of 10^. The coherent component of the radiation, i.e., the 
sharply focused subbeams that decay as -Rp^^^, are emitted at the frequency /xo; The 
spherically decaying, incoherent component of the radiation arising from the polarization 
current described in Eq. ([T]), on the other hand, contains frequencies that are higher than 
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/icu by a factor of order (yt/ujY ■'■^ pulsars, (yt/ujY ~ 10^^ when Vt/uo is comparable to 
/i, z.e., when the frequency moj that characterizes the spatial fluctuations of the emitting 
plasma is of the order of, or smaller than, its modulation frequency Vt. Hence, not only the 
power-law indices of the coherent component, but the unusually broad spectral distribution 
of the incoherent component of this radiation, too, is consistent with the observational data 
from certain pulsars. The pulsed emission from the Crab pulsar, for example, extends over 
53 octaves of the electromagnetic spectrum from radio waves to 7-rays [l5j. 

We note, finally, that neither the asymptotic expansion presented here nor that which 
was obtained in Refs. [3, B| are uniform with respect to the parameter 9-p. The present 
approximation receives contributions only from the volume elements in the vicinity of the 
single source point f = 1, z = ip at which the cusp curve A = of the bifurcation surface 
touches the light cylinder (see Figs. 3 and 4 of Ref. [?!). This is in sharp contrast to 
the asymptotic expansion of the field outside the equatorial plane, where the leading term 
receives contributions from a filamentary locus of source elements, i.e., from the intersection 
of the cusp curve of the bifurcation surface with the entire volume of the source. Comparison 

n 

of Eq. ( 1271) with its counterpart, Eq. (57) of Ref. [7|, shows that the (smooth) transition 
from the nonequatorial to the equatorial regime occurs across 9-p ^ arccos(yU^^/^). However, 
the derivation of a uniform asymptotic approximation to integral 1 that would determine 
the field in the transition region is a challenging mathematical problem that remains open. 
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